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Abstract 

A liquid Argon TPC (ICARUS-like) has the ability to detect neutrino bursts from 
type-II supernova collapses, via three processes: elastic scattering by electrons 
from all neutrino species, and Uf. charged current absorption on Ar with produc- 
tion of excited K and charged current absorption on Ar with production of 
excited CI. In this paper we have used new calculations in the random phase 
approximation (RPA) for the absorption processes. The anti-neutrino reaction 
Vf. '^^ Ar absorption has been considered for the first time. For definiteness, we 
compute rates for the 3 kton ICARUS detector: 244 events are expected from a 
supernova at a distance of 10 kpc, without considering oscillation effects. The 
impact of oscillations in the neutrino rates will be discussed in a separate paper. 
We also discuss the ability to determine the direction of the supernova. 
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1 Introduction 



When a massive star (eight or more times the mass of the Sun) runs out of nuclear 
fuel, it develops a layered structure, which heavier elements lay in deeper layers, the 
innermost element being iron. This configuration becomes eventually unstable and 
leads to a catastrophic collapse of the star core in a time window of a few thousandths 
of a second, since the star cannot resist the pressure of its internal gravitational force. 
This collapse leads to an explosion that is known as type II supernova [T] . The expected 
rate for such an event to occur in our Galaxy is one every thirty years j2]. 

The binding energy released during this process is enormous: 10^^ ergs. The high 
density prevents the emission of photons, therefore 99% of the energy comes out in 
the form of neutrinos of all types. A small fraction of them, about 1%, are produced 
during the neutronization process through the reaction e~ + p ^ lye + n, while the rest 
are uu pairs from later cooling reactions. 

On their way out, neutrinos interact inside the stellar core which is an extremely 
dense and neutron-rich medium. The different scattering reactions, underwent by ev- 
ery neutrino species, result in different energy distributions for the emitted neutrinos. 
Muon and tau neutrinos interact only through weak neutral currents, thus they decou- 
ple deep inside the core, where a higher temperature exists, and therefore are emitted 
with higher energies compared to z/g's and UeS. In addition, z/g's suffer fewer charged 
current interactions than i/gS, so they decouple in hotter regions, thus having an energy 
spectrum that is harder than that of z/g's. 

The neutrino signal from a supernova rises first steeply and then falls exponentially 
with time (see The experimental detection of a neutrino signal from a supernova 
is therefore expected to occur in a time window of about 10 seconds. This implies that 
detector thresholds can, in principle, be set quite low since backgrounds are highly 
suppressed given the narrow time window in which the burst takes place. 

The detection of neutrinos emitted by SN1987A in Kamiokande |1] and 1MB [S] 
began the era of neutrino astronomy. Nowadays, several running or planned experi- 
ments (SuperKamiokande, SNO, LVD, KamLAND, Borexino...) have the capability to 
unambiguously detect supernova neutrinos [Ul |7] . 

As discussed in jHIH], a liquid Argon TPC like in ICARUS has the ability to detect 
neutrino bursts from supernova collapses. Preliminary estimates of the expected rates 
yielded 76 events (44 from absorption and 32 from electron scattering assuming a 
5 MeV electron detection threshold) for a supernova at 10 kpc. For the absorption 
cross section, a threshold of 11 MeV was taken into account^. The energy threshold 
for this Fermi transition is however 5.885 MeV [9]. In addition, the Gamow- Teller 
transitions were not included in the total absorption rate. 

In this paper, we improve the estimates for supernova neutrino rates taking into 

^in ref. ^ a threshold Q — 11.5 is quoted in units of rrie. 
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account latest developments: we use new calculations in the random phase approx- 
imation (RPA) for the absorption processes JU] and consider all the processes that 
contribute to the total number of expected neutrino events, including for the first time 
the absorption reaction v^. ^^Ar. We discuss the ability to determine the direction of 
the supernova. 



2 Supernova neutrino signatures 

Models of type II supernovae predict that neutrinos are emitted with a thermal spec- 
trum, with a temperature hierarchy among neutrino flavours: Ty^ < Tp^ < Ty^ 
The neutrino energy spectra can be described by a Fermi-Dirac distribution |llj : 

dN C El 



where C = 0.55; T is the temperature (MeV); E^, the neutrino energy (MeV); rj the 
chemical potential and the number of expected neutrinos of a given species. We 
use r] = ^T]. The following values for the different neutrino temperatures are used, 
giving the following average energies: 
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According to [12], pairs of neutrino flavours should be produced with equal luminos- 
ity. The total number of neutrinos of a given flavour is obtained dividing the binding 
energy by six (corresponding to all + 1> species) and by the average energy associated 
to each species. As a result, will be more copiously produced than u^^s since 
their mean energy is lower. Figure ^ shows the expected neutrino spectra for a binding 
energy of 3 x 10^^ ergs. We expect 2.8 x 10^'' z/g, 1.9 x 10^'' z/g and 5.0 x 10^'^ i/^^r + i^fi,T 
for a total of 9.7 x 10^^ neutrinos. 

In argon two processes contribute to the total rate: 

• Elastic scattering: z^x + e~ ^ z/^ + e~ (x = e, /i, r) sensitive to all neutrino 
species. 

• Absorption: Ar e" K*, Ar e+ CI* 

We consider two configurations of the detector with masses of 1.2 and 3 kton: the 
number of targets is accordingly: 

1.2 kton 3 kton 

Elastic scattering 3.2 x 10^^ e~ 8.0 x 10^^ 

Absorption 1.8 x 10^^ Ar nuclei 4.5 x 10^^ Ar nuclei 
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2.1 Elastic events 



The elastic neutrino scattering off electrons has a total cross section that increases 
linearly with energy: 

a{iyee- Uee') = 9.20 x IQ-^^ E^^{MeV) cin 

a{vee- Vee ) = 3.83 x IQ-^^ Ep^{MeV) cin , . 

a{iy^^,e~ ^ u^^rC-) = 1.57 x IQ-^s^,^ ^ (MeV) cm' 

(T(/>^,,e- z/^,,e-) = 1.29 X 10-45 Ep^ ^ (Me V) cm' 



All neutrino species contribute to elastic scattering. The experimental signature 
consists in a single recoil electron. Since, the direction of this electron is highly cor- 
related to the incoming u direction, these events have the potentiality of precisely 
determining the location of the supernova source. 

Table^shows the expected rates for this reaction. To compute the number of events 
we fold the total cross section as a function of energy with the appropriate Fermi-Dirac 
distribution. Our rates are calculated integrating over all electron recoil energies. Since 
the neutrino burst occurs in a time window of about 10 seconds, we estimate that the 
background expected due to natural radioactivity is negligible and therefore we do not 
apply any threshold for electron detection. 

Figure |21 shows the expected elastic event rates, as a function of the incoming 
neutrino energy, for the case of a 3 kton detector and a supernova occurring at 10 kpc. 
The largest contribution (~ 50%) to elastic events comes from z/g, since they have both 
the larger flux and cross section. For a 3 (1.2) kton detector, we expect a total of 41 
(16) elastic events. 
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40 Ar 


3.5 


11 


75 
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40 Ar 


5 


16 


6 


15 


Total 








97 


244 



Table 1: Expected neutrino rates for a supernova at a distance of 10 kpc, releasing an 
energy of 3 x 10^^ ergs (no threshold on the electron energy has been applied). 
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2.1.1 Reconstruction of the incoming neutrino direction 

It is possible to use the angular distribution of elastic scattering neutrino events to 
identify the direction of a nearby supernova. Ordinarily, this would be unnecessary due 
to the optical emissions. However, if one is interested in studying a galactic supernova's 
first hours or days or if the supernova is obscured in the optical by the dust of the 
galactic disk, the early warning and crude angular position provided by the neutrino 
data would be uniquely useful. 

One method to achieve this is to use the fact that the electron produced via an 
elastic scattering is strongly forward-peaked and it will point back to the supernova 
direction. Depending on the type of the incident neutrino and its energy, one can 
expect the distributions of the secondary electrons to be confined to a certain effective 
scattering cone. 

In order to estimate our angular efficiency, supernova neutrino events have been 
fully simulated using the FLUKA Monte Carlo package [12]. A detailed description of 
the ICARUS T600 geometry is included in the simulation ^3]. The detector internal 
volume consisted of two half-modules filled with liquid argon. A fiducial volume of 3.2 
X 3.0 X 18.0 m^ has been considered for every half-module, surrounded by a dewar of 
stainless steel of 20 tons. An external container of aluminum and honeycomb acts as 
thermal insulator. Finally, a 70 cm thick layer of polyethylene tubes filled with boric 
acid has been simulated as a shield against neutrons. 

Elastic events are characterized by the production of a primary electron track in the 
sensitive volume in a direction correlated with the incoming neutrino direction. The 
electron undergoes multiple scattering and Bremsstrahlung which leads to a deviation 
from the original direction, also affecting the angular reconstruction accuracy. 

The signal deposited by the track in the detector is discretized in elementary cells, 
which correspond to the readout granularity of the detector. The detector volume is 
supposed to be divided in cells of 3 x 3 x 3 mm'^. Then, the energy associated to the 
track is computed from the charge deposited in every elementary cell. 

The reconstruction of the electron momentum is performed by using the position of 
every cell with respect to the primary vertex weighted with the corresponding charge 
deposited in the cell. The algorithm can be optimized by varying the track length 
(and the consequent number of cells) considered in the reconstruction of the electron 
direction to avoid the multiple scattering effect. 

Two different optimization methods have been tested in terms of two variables: the 
percentage of cells belonging to the electron track to be used in the reconstruction and 
the scattering angle of the electron. 

In the first method we studied the variation of the angular resolution as a function 
of the percentage of cells taking into account in the momentum reconstruction. The 
angular resolution was computed as the RMS of the distribution ()rec_QMC ^ ^.j^g 
ference between the reconstructed and simulated angles between the electron and the 
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parent neutrino directions. Figure |3 shows this variation for different electron energy 
ranges. The percentage of cells used in the electron direction reconstruction giving the 
best angular resolution decreases when the electron energy increases. 

The reconstructed and simulated angles for different electron energy ranges are 
plotted separately in figure El The value of the angle between the electron and the 
neutrino is dominated by the kinematics of the events and not by the reconstruction 
algorithm. 

Figure El shows the expected angular resolution as a function of the electron energy 
for different percentage of cells. For low electron energies (< 2.5 MeV), almost all cells 
are needed to have an accurate reconstruction (90%). For energies between 2.5 and 16 
MeV, the optimal resolution is obtained using the 50% of the cells and for electrons of 
more than 16 MeV only the 30% of the cells should be used. 

Figure [7\ shows the angular resolution obtained with the best reconstruction algo- 
rithm which uses the optimal percentage of cells for each electron energy range. The 
resolution decreases strongly at low energies, being below 16 degrees for events with 
electron energies of more than 5 MeV. For energies greater than 15 MeV the resolution 
is less than 8 degrees. 

The second method is based on the dependence between the electron momentum 
and the scattering angle. The projected scattering angle can be obtained as 15]: 



where x is the track length, p and Pc are the momentum and velocity of the particle, 
and Xq = 14 cm is the radiation length in LAr. Therefore the particle momentum, 
neglecting ionization losses, is proportional to the inverse of the average value of the 
scattering angles. 

We select the track length for which the scattering angle is smaller that a cer- 
tain threshold. This way, we reduce the effect of multiple scattering on the angular 
resolution. 

Figure |H1 shows the angular resolution as a function of the electron energy for dif- 
ferent values of the scattering angle threshold. The same behaviour of the angular res- 
olution as a function of the electron energy is observed with this method but slightly 
worse results are obtained. Events with electrons of more than 5 MeV have an an- 
gular resolution below 18 degrees. For more than 15 MeV, the resolution reached is 
essentially the same with both methods. 

Therefore, the first algorithm based in the percentage of cells has been adopted for 
the reconstruction of the scattered lepton direction. 

Figure El shows the fraction of events as a function of the cone angle within which 
the reconstructed electron direction in space is contained. The cone axis is defined by 




(3) 
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the parent neutrino direction. The different curves correspond to different thresholds 
on the electron kinetic energy (Tiep)- For a 25° cone around the neutrino direction and 
electron energies greater than 2 MeV, the angular efficiency is about 51%. It increases 
as a function of the threshold on the electron energy, being 66% for electrons of more 
than 5 MeV. 

In order to determine the supernova explosion location, the angular distribution 
of the elastic scattering neutrino events was used. Figure ^1 shows the distribution 
of the angle between the incoming neutrino and the outcoming electron for different 
electron energy thresholds. The uncertainty in the direction can be computed from 
the RMS of the distributions shown in figure ^2 The sign of the angle refers to 
the projection of the electron momentum into an axis perpendicular to the neutrino 
direction. The contamination coming from absorption events reconstructed as elastic 
events is indicated in the plot. It is largely reduced as soon as we apply a cut on the 
lepton energy. 

Therefore, the resulting uncertainty in the direction measurement of a supernova at 
a distance of 10 kpc is ~ 70.6° /^/N for Ti^p > MeV, ~ 38.97v^ for T,ep > 2.5 MeV 
and ~ 28.8°/a//V for T/ep > 5 MeV, being N the expected number of elastic supernova 
events. The corresponding uncertainties obtained with a 1.2 kton and 3 kton detector 
are given in table El together with the expected elastic rates for different lepton energy 
thresholds. 



T 

-■- lep 


Expected elastic events 


Expected uncertainty 


(MeV) 


1.2 ktons 


3 ktons 


1.2 ktons 


3 ktons 


> 


16 


40 


17.7° 


11.2° 


> 2.5 


12 


30 


11.2° 


7.1° 


> 5 


9 


23 


9.6° 


6.0° 



Table 2: Expected uncertainty in the direction measurement of a supernova at a dis- 
tance of 10 kpc for a 1.2 kton and 3 kton detector for different electron energy thresh- 
olds. 



2.2 Absorption events 

The h'e'^^Ar absorption rate is expected to proceed through two main channels: a 
superallowed Fermi transition to the 4.38 MeV excited isobaric analog ^°K* state; 
Gamow- Teller transitions to several excited states. These two processes can be 
distinguished by the energy and multiplicity of the 7 rays emitted in the de-excitation 
and by the energy spectrum of the primary electron. In addition one can consider other 
processes which become important energies above 20 MeV, where the final state is not 
a bound state anymore and where emission of particles from the nucleus is possible. 
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The cross section for the Fermi capture is given by jH]: 



o = 1.702 X IQ-^^^^ E^-mlF{Ee) cm^ (4) 

where F{Ee) is the Fermi function, Ee = — Q + me is the prompt electron energy, 
Q is the energy threshold [Q = 5.885 MeV) and nie is the electron mass. 

In reference jTHj shell model calculations of Ug absorption cross sections for '^^Ar 
are presented by Ormand et al. They show that GT transitions lead to a significant 
enhancement of the absorption rate over that expected from the Fermi transition. 

In recent studies |10] the z/e absorption cross section has been computed by random 
phase approximation (RPA) for neutrino energies up to 100 MeV. These calculations 
include all the multipoles up to J=6 and both parities, whereas only the 1"*" and O"*" 
allowed transitions were considered in Ref. |l6j. As shown in figure El both results 
are similar up to ~ 20 MeV. At higher energies relevant for supernova neutrinos, the 
cross section including all the multipoles increases faster, showing the importance of 
incuding higher multipoles in order to get a better estimate of the event rate. 

The absorption cross section is much larger than that for v^. elastic scattering, hence 
this process significantly enhances the experimental sensitivity (see table P). In a 3 
kton detector, we expect 188 Ve absorption events from the supernova at 10 kpc. 

The absorption process +^'^ ~^ +^'' CI* is included for the first time, giving 
15 events for a 3 kton detector. The corresponding cross section computed by RPA [TU] 
is also plotted in figure El In this case the l"*" and O"*" allowed transitions are suppressed 
by nuclear structure reasons and the cross section is dominated by the other multipoles. 

The total expected rate for a 3 kton (1.2 kton) detector is about 244 (97) events. 
This statistics would allow a detailed study of neutron star and even black hole for- 
mation and would largely contribute to get an improved knowledge of some neutrino 
intrinsic properties like, for example, the mass, electric charge and magnetic moment. 

Figure ^1 show the total rate for a 3 kton detector as a function of supernova 
distance. It also displays the expected rates for the different neutrino reactions. For 
a gravitational stellar collapse occurring in the Large Magellanic Cloud (distance ~ 
60 kpc), we expect to collect, in a 3 kton detector, about 7 events as a result of the 
neutrino burst. 

3 Conclusions 

A liquid Argon TPC like ICARUS is sensitive to supernova neutrinos. A typical gravi- 
tational stellar collapse, occurring at 10 kpc, will produce a neutrino burst of about 97 
(244) events in a 1.2 (3) kton detector. The largest contribution comes from the ab- 
sorption process of Pe neutrinos. Elastic scattering off electrons for all neutrino species 
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contributes to 15% of the total rate. The Vf, absorption reaction has been studied for 
the first time, giving 15 events for a 3 kton detector. 

The direction of the lepton momentum for elastic scattering neutrino events can 
be measured with good accuracy. A reconstruction algorithm has been optimized by 
varying the track length to be considered in the reconstruction. The contamination 
from absorption events in the elastic channel can be reduced cutting on the lepton 
energy. The expected angular resolution is below 16 degrees for lepton energies greater 
than 5 MeV. The resulting uncertainty in the direction measurement of a supernova 
at 10 kpc is 17.7° (11.2°) with a 1.2 (3) kton detector, reaching 9.6° (6.0°) for lepton 
energies greater than 5 MeV. 

Neutrino oscillations will change significantly the expected number of events. Mat- 
ter effects in the SN and vacuum oscillations will modify the neutrino fluxes. As a result 
of mixing, the original v^. spectrum will be shifted to higher energies. This will increase 
dramatically the rate of Ug^^Ar absorption events due to the strong dependence of the 
cross section with the neutrino energy. Because of its sensitivity to neutrinos, a 
liquid Argon TPC can also provide unique information about the early breakout pulse 
that might be observed from the next galactic supernova. Analyses considering these 
features will be included in a following paper [T7j. 
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Supernova neutrino energy spectra 
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Figure 1: Supernova neutrino energy spectra for a type II supernova releasing 3 x 10 
ergs. 
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Supernova neutrino rates in 3 kton ICARUS 
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Figure 2: Elastic (top) and absorption (bottom) neutrino event energy distribution for 
a type II supernova at a distance of 10 kpc. 
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Figure 3: Cross sections as a function of neutrino energy for z/g and Pe absorption 
reactions. We compare the results from shell model calculations fB] considering only 
Fermi and Gamow-Teller transitions with those from RPA calculations including all 
the multipoles [TU] . 
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Figure 4: Supernova neutrino angular resolution as a function of the percentage of cells 
used in the momentum reconstruction for different lepton energy ranges. 
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Figure 5: Reconstructed {6^'^'^) and generated [9^^^) angles between the leading electron 
and the incoming neutrino as a function of the percentage of cells used in the momentum 
reconstruction for different lepton energy ranges. 
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Angular resolution 
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Figure 6: Supernova neutrino angular resolution as a function of the lepton energy for 
different percentage of cells used in the momentum reconstruction. 
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Figure 7: Supernova neutrino angular resolution as a function of the lepton energy. 
The optimal algorithm for the momentum reconstruction is used. 
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Figure 8: Supernova neutrino angular resolution as a function of the lepton energy for 
different values of the scattering angle threshold. 
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Supernova neutrinos 
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Figure 9: Fraction of events as a function of the angle between the scattered electron 
and the parent neutrino direction. The different curves correspond to different lepton 
energy thresholds. 
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Figure 10: Angle between the incoming neutrino and the outcoming electron for elastic 
events with electron energies greater than 0, 2.5 and 5 MeV, respectively. 
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Mean 0.1569 
RMS 28.82 
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Figure 11: Angle between the incoming neutrino and the outcoming electron for elastic 
events with electron energies greater than 0, 2.5 and 5 MeV, respectively. The sign of 
the angle refers to the projection of the electron momentum into an axis perpendicular 
to the neutrino direction. 
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Figure 12: Expected neutrino rates for elastic and absorption reactions as a function 
of the supernova distance. We consider a 3 kton ICARUS-like detector. 
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Supernova neutrino energy spectra 
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